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Abstract 

The most general one dimensional reaction-diffusion model with nearest- 
neighbor interactions solvable through the empty interval method, and 
without any restriction on the particle-generation from two adjacent empty 
sites is studied. It is shown that turning on the reactions which generate 
particles from two adjacent empty sites, results in a gap in the spectrum 
of the evolution operator (or equivalently a finite relaxation time). 
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1 Introduction 



Reaction-diffusion systems have been studied using various methods, includ- 
ing analytical tcKihniqucs, approximation methods, and simulation. The proper 
approximation methods are generally different in different dimensions, as for 
example the mean field techniques, working good for high dimensions, generally 
do not give correct results for low dimcmsicmal systems. A large fraction of an- 
alytical studies, belong to low-dimensional (specially one-dimensional) systems, 
where solving low-dimensional systems should in principle be easier. [1-13]. 

In this context, the term solvability (or integrability) is used in different 
senses. In [14-16], integrability means that the A^-particle conditional probabil- 
ities' S-matrix is factorized into a product of 2-particle S-matrices. In [17 27], 
solvability means closedness of the evolution equation of the empty intervals (or 
their generalization). In [28-30], solvability means that the evolution equation 
of n-point functions contains only n- or less- point functions. 

Among the important aspects of reaction-diffusion systems, are the station- 
ary state of the system (or one of the quantities describing the system) and the 
relaxation behavior of the system towards this configuration. In the thermody- 
namic limit (when the size of the system tends to infinity) these behaviors may 
show discontinuity in terms of the control parameters of the system. In [31-34] 
(for example); such behaviors are studied. 

The empty interval method (EIM) has been used to analyze the one dimen- 
sional dynamics of diffusion-limited coalescence [17-20]. Using this method, 
the probability that n consecutive sites are empty has been calculated. This 
method has been used to study a reaction-diffusion process with three-site in- 
teractions [22] . EIM has been also generalized to study the kinetics of the q'-state 
one-dimensional Potts model in the zero-temperature limit [21]. 

In this article, we are going to study all the one dimensional reaction-diffusion 
models with nearest neighbor interactions which can be exactly solved by EIM. 
It is worth noting that ben-Avraham et al. have studied one-dimensional 
diffusion-limited processes through EIM [17-20]. In their study, some of the 
reaction rates have been taken infinite, and they have worked out the models on 
continuum. For the cases of finite reaction-rates, some approximate solutions 
have been obtained. 

We study models with finite reaction rates, obtain conditions for the system 
to be solvable via EIM, and then solve the equations of EIM. In [23], general con- 
ditions were obtained for a single-species reaction-diffusion system with nearest 
neighbor interactions, to be solvable through the empty-interval method. Solv- 
ability means that evolution equation for Ep (the probability that n consecutive 
sites be empty) is closed. It turned out there, that certain relations between the 
reaction rates are needed, so that the system is solvable via EIM. The evolution 
equation of En is a recm-sive equation in terms of n, and that this equation 
is linear. It was shown that if certain reactions are absent, namely reactions 
that produce particles in two adjacent empty sites, the coefficients of the empty 
intervals in the evolution equation of the empty intervals are n-independent, 
which makes them be solved more easily. The criteria for solvability, and the 
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solution of the empty-interval equation were generalized to cases of multi-species 
systems and multi-site interactions in [24,26,27]. 

Here we want to study the case dropped from the study in [23], namely when 
there arc interactions producing particles from two adjacent empty sites. Doing 
so, we are considering the most general one dimensional reaction-diffusion model 
with nearest-neighbor interactions which can be solved exactly through EIM. 

The scheme of the paper is as follows. In section 2, the most general one 
dimensional reaction-diffusion model with nearest-neighbor interactions which 
can be solved exactly through EIM is introduced. In the same section the evo- 
lution equation of the empty intervals is obtained for a lattice. Then, using a 
limiting procedure a similar equation is obtained for the continuum. In section 
3 the stationary solution to this equation is obtained. In section four the re- 
laxation of the system towards its stationary state is investigated. Section 5 is 
devoted to the concluding remarks. 

2 Models solvable through the empty interval 
method 

To introduce the notation, let us briefly review the criteria that a single-species 
nearest-neighbor-interaction reaction-diffusion system be solvable through the 
empty- interval method (EIM). Consider a one-dimensional lattice. It was shown 
in [23], that the most general interactions for a single-species model in a one- 
dimensional lattice with nearest-neighbor interactions are 



in order that the system be solvable through the EIM. Here an empty (occu- 
pied) site is denoted by o (•), and r^'s and r are reaction rates. Denoting the 
probability of finding n consecutive empty sites by 




(1) 



and 



oo ^ anything, 



(2) 



n 




o o • • • o) =: E, 



(3) 



it was then shown that 



dt 



= {r2 + r3){En-l + En+l - 2 

- {n-l)rEn, n > 1, 



2En) - in + r4){En - En+i) 



(4) 



dEijt) 
dt 



=(r2 + r3)(l + E2- 2Et) - (n + ri){Er - E2) 



(5) 



dEL+i{t) 
dt 



= - LtEl+i, 



(6) 
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where the length of the lattice has been assumed to be L + 1 . It is seen that 
the (O takes a form similar to provided one defines 

Eo{t) 1. (7) 

In [23], 10} to © were actually obtained for the case r = 0. Q to iQ, and of 
course the initial values of i?„'s, are a complete set of equations to obtain En{t). 
One can absorb the rate (r2 + ^3) in the definition of time, and rewrite Q and 
(0 as 

7^-^ —{En-l + En+l — 2En) ~ b (£'„ — En+l) 

at 

-{n-l)cEn, 0<n<L + l, (8) 



where 



b — , 

r2 + rz 
r 

c :=- 



f2 + rs 



(9) 



The aim is to solve (jHl along with ((HJ and 

The continuous-space form of the above equations is 

^-^ + b^~c.E, Q<.<X, (10) 
ot ox'' ox 

with the boundary conditions 

E{x^O,t)^l, (11) 
?E>£^ = -r.XEi. = X.r, (12) 

The procedure to write these equations in continuous space is to define 

X := n A, 

c 

E{x,i) -.^ En{t). (13) 

One then expands the right-hand sides of and JHJ in terms of A, sends A to 
zero, and substitutes the quantities with tilde with the corresponding quantities 
without tilde. 

Using the new variable £ defined through 

£{x,t):^E{x,t)exp(^Y (14) 
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one can rewrite (|1U|) to (|12|l as 



^-a^:^- 4+-)^. o<-<x, 



S{x = 0,i) = 1, 
dS{x = X,t) 



cX£{x = X,t). 



(15) 
(16) 
(17) 



3 The stationary solution 



Denote the stationary solution to (|10|l and (|12() by . It is seen that £^ is a 
hnear combination of the Airy functions, so, 



i?P(x-)=exp(-^) ^aAi 



c I ex - 



/3Bi 



c I C2; + 



where a and /? are two constants satisfying 

'c-2/3 52 



(18) 



fAi 



4 

-/3Bi 



PBi{^-^ I -1, 



,-2/3 



ca: + — 



= 0. 



(19) 



This solution is simplified for X —^ 00 (the thermodynamic limit). As Bi(?/) 
behaves like the exponential of a y'^^^ for y ^ 00, it is seen that in the thermo- 
dynamic limit P is zero. So, 



1 /_bx ^. 



c I ex H 



a: ^ 00. 

(20) 

It is seen that there is a unique stationary solution. 

The above argument is valid for c 7^ 0. If c = 0, then E{x — X,t), and hence 
£{x — X,t), is i-independent, and one has 



£;^(a;) = 7 + (1 - 7) exp(-6x), c = 0, 



(21) 



where 7 is an arbitrary constant between zero and one. It is seen that in this 
case the stationary solution is not unique. It is also noteworthy that (|21|l is not 
the limit of H18() or (|20|l as c tends to zero. 



4 Relaxation towards the stationary solution 

Defining 

F{x, t) := E{x, t) - E^{x), (22) 
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it is seen that the evolution equation for F is the same as that of E, except 
for the fact that the boundary conditions for F are homogeneous. To calculate 
F{x,t), one seeks the eigenvalues and eigenvectors of the evolution operator: 



^Je{x) ^-7-^ +0- ca;/,, 

dx^ da; 

/e(0) =0, 

The solution to this is 
fe{x) =exp ( 



0<x<X, 



(23) 

(24) 
(25) 



X < a Ai 



,-2/3 



cx + e + 



/3Bi 



,-2/3 



CX + e + 



(26) 



where where a and [3 are two constants satisfying 

&2M 



,-2/3 



e + 



/3Bi 



-2/3 



!Ai 



,-2/3 



cX 



4 



/3Bi 



-2/3 cX 



62 
4 
62 
4 



-0, 
= 0. 



(27) 



It is seen that the above equations for a and (3 have nonzero solutions, only 
for certain discrete values of e. This means that the spectrum of the evolution 
operator is discrete, and there is a gap between zero and the largest nonzero 
eigenvalue of the evolution operator. Here too, the solution is simplified if one 
considers the thermodynamic limit. In this case. 



f^{x) = exp 
where e is among e„'s: 



hx 
T 



Ai 



-2/3 



CX 



62 
4 



X 



C ^ Zfi 



4 ' 



and Zn 's are the zeros of the Airy function: 

Ai(z„) = 0. 



(28) 

(29) 
(30) 



It is seen that if one tends c to zero and X to infinity, the spectrum of the 
evolution operator tends to (— oo,— 62/4). However, if one puts c = and 
X = oo, and then solves the eigenvector equation, another result is obtained. 
In this case, (^5)1 to become 



eUx) 



d2/. 



dx^ 

MO) =0, 

lim /,(x) =0. 



6^, 0<x, 
ax 



(31) 

(32) 
(33) 
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The solution to these is 



fe{x) = sinh ( y e + 2; I exp 





) 



(34) 



and the only condition for e is that e must be negative. That is, the spectrum 
of the evolution operator is (—00, 0). 

5 Concluding remarks 

The most general one-dimensional single-species exclusion model was consid- 
ered, for which the evolution of the empty-intervals is closed. The effect of 
particle creation in two empty adjacent sites was specially investigated. The 
stationary solution was obtained and the relaxation towards this stationary so- 
lution was studied. It was shown that if the rate of particle creation in adjacent 
empty sites is nonzero, then the spectrum of the evolution operator of the empty 
intervals is discrete. If this rate is zero and the system is infinite, then the spec- 
trum is continuous. However, the spectrum depends on whether one finds the 
spectrum for the finite system and then tends the size of the system to infinity, 
or the spectrum is directly calculated for the infinite system. In the former case, 
the largest eigenvalue of the evolution operator is negative (there is a gap in the 
spectrum) and the results of [23] are recovered. This means the the relaxation 
of the system towards its steady state is exponential, in other words, the system 
has a finite relaxation time. In the latter case, there is no gap in the spectrum 
and the spectrum extends to zero. So in this case the relaxation of the system 
towards its steady state is not exponential, in other words, the relaxation time 
of the system is infinite. This is an example of a system for which the limit 
of the spectrum as the size of the system tends to infinity is different from the 
spectrum of the infinite system [1]. 
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